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Abstract 

We investigate the Fermi Liquid theory of high density QCD. Using the renormal¬ 
ization group (RG), we determine the behavior of effective fermion interactions near 
the Fermi surface. At sufficiently high densities the matching between the Fermi Liquid 
theory and QCD can be accomplished in perturbation theory and a reliable calculation 
of Cooper pair formation performed, modulo the existence of a magnetic screening 
mass for the gluon. The presence of a chemical potential leads to different RG flows 
for sub-components of the effective interactions which would ordinarily be linked by 
Lorentz invariance. We also study instanton-induced operators and find that near the 
Fermi surface they are likely to be subdominant relative to operators induced by per¬ 
turbative gluon exchange. We discuss the implications of our results for the phase 
structure of QCD at finite baryon density and temperature. 
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1 Introduction 


In this paper we investigate the dynamics of quark excitations near the Fermi surface of 
high density QCD. Asymptotic freedom, as well as explicit calculations |l|, suggest that 
the high density limit is controlled by a small coupling: evaluated at the length scale 

of typical quark separation. However, it is well known that even very weak interactions 
can lead to interesting behavior in the presence of a Fermi surface P|. In particular, pair 
formation leading to a mass gap, superfluidity and superconductivity seem generic whenever 
an attractive interaction, no matter how weak, is present. One way to understand this is 
through the Fermi Liquid theory originally developed by Landau, and recently clarihed using 
the renormalization group (RG) We review the RG method in section 2 below within 
the context of relativistic theories. Our paper can be regarded as a first exploration of this 
method applied to ultra-relativistic, non-Abelian gauge theories. Here we simply note that 
Fermi Liquid theory can be regarded as an effective held theory of fermionic excitations, in 
which all interactions turn out to be irrelevant with the exception of four fermion vertices 
with the kinematic structure associated with Gooper pairing (i.e. - two quarks scattering 
with equal and opposite spatial momenta). The renormalization group hows of these special 
operators therefore determine the characteristics of the model. 

At sufficiently high density we expect that the matching between the fermionic ehective 
theory and QGD can be accomplished in perturbation theory (we discuss the details of this 
matching in section 3). The four fermion operators which are induced at leading order are 
simply those from one gluon exchange with an IR cut oh provided by the gluonic screening 
masses. The gap equations resulting from one gluon exchange were studied in some detail 
by Baffin and Love Q]. As we will discuss in section 4, their results are qualitatively similar 
to ours in the high density regime, with the attractive 3 interaction leading to 0^ or 0“ 
diquark condensation and color superconductivity. The RG approach provides an intuitive 
understanding of this result since it is the attractive 3 interaction between two left (right) 
handed quarks that hows hrst to a Landau pole as we approach the Fermi surface. 

At somewhat lower densities, one might expect that non-perturbative interactions, such 
as those mediated by instantons, play a role in the Fermi Liquid theory. Recent attention 0 
has focused on the ehects of instanton-inspired fermionic models of QGD, although most of 
the analysis has been in the mean-held approximation, with quantum corrections neglected. 
In section 5 we study the combined RG how of instanton and one-gluon induced vertices. 
We hnd that the one-gluon vertices grow much faster as the Fermi surface is approached. 
Therefore, unless the instanton vertices are larger than the one-gluon vertices at the matching 
scale (which seems implausible), the dynamics is probably controlled by one-gluon ehects. 
This suggests that the high density phase identihed above will persist to the (assumed) 
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boundary with the low density phase which exhibits chiral symmetry breaking. 

Our paper is organized as follow. In section 2 we discuss the RG description of Fermi Liq¬ 
uid theory for relativistic systems. In section 3 we describe the matching conditions between 
QCD and our effective theory. In section 4 we derive and solve the RG equations governing 
the flow of gluonic four fermion operators, obtaining an expression for the condensate size 
as a function of quark density. In section 5 we study the combined RG equations for the 
instanton and gluonic operators, and hnd that the gluonic operators are dominant near the 
Fermi surface. In section 6 we discuss our results and their implications for QGD at high 
densities and temperatures. 


2 Fermi Liquids from the Renormalization Group 

We wish to study QGD at finite density, which is described by the following lagrangian with 
chemical potential /i 

£ = + A(ip+ n)A ■ (i) 

We make a guess as to the form of the effective theory close to the Fermi surface; the 
obvious guess based on the dynamics of non-relativistic systems is that the theory is one 
of weakly interacting fermions: these are the dressed “quasi-particles” of solid state physics 
language. We will henceforth refer to these effective degrees of freedom as “quarks”, with 
the understanding that they may in fact be related to the bare quarks in a complicated way. 
Rather than treating the gluons as propagating degrees of freedom we will integrate them 
out leaving a potentially infinite sum over local, higher dimension fermion operators. The 
locality of these operators requires that gluons be screened at long-distances, presumably 
by effective electric and magnetic mass terms induced by the medium. We will address the 
matching between the high scale QGD dynamics and the four Fermi theory in section 3 but 
for the moment we assume the matching can be performed and proceed. 

The Fermi surface in (^) naturally breaks the 0(3,1) invariance of space-time to 0(3) 
and furthermore picks out momenta of order fi. It is therefore natural to study the theory 
as we approach the Fermi surface in a Wilsonian sense. We parameterize four momenta in 
the following fashion 

= {Po,P) = {K.k + l) ( 2 ) 

where k lies on the Fermi surface and I is perpendicular to it. We study the Wilsonian 
effective theory of modes near the Fermi surface, with energy and momenta 

|A:o|,|/l < A , A^O . (3) 
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While this type of RG scaling is somewhat unfamiliar, it actually corresponds to thinning 
out fermionic degrees of freedom according to their eigenvalues under the operator + /lyo- 
It is easy to see that eigenspinors of this operator with eigenvalues Xn ■ |An| < A correspond 
to states satisfying Consider an eigenspinor of the form 

= u{E,p) , (4) 

where u{E,p) satisfies the usual momentum-space Dirac equation with E = \p\ = \k + 1\ = 
/i ± 0{A), and po = 0{A). Then by direct substitution we see that xjjp satishes 

{i(^ + /i7o) V’fc = C>(A) . (5) 

Thus, the RG flow towards the Fermi surface just corresponds to taking the cutoff on eigen¬ 
values of + p'jo to zero. 

Which interactions are relevant operators in this limit? For our guess to make sense 
the kinetic term for the fermions must be invariant when we scale energies and momenta, 
ko —> sko (or, t —>• t/s), and I —>• si, with s < 1. We must be careful though to satisfy all the 
global symmetries of the theory. In particular, there is a spurious symmetry of (|I]) in which 
we treat p as the temporal component of a fictitious U{1)b gauge boson. In other words, 
the combined transformations 'ip —>■ e*®* pj and p ^ p — 6 leave the lagrangian (^) invariant. 
From this, we deduce that derivatives acting on ip may only enter the effective theory as 
i(p -|- /i 7 o. This requires the kinetic term of our effective theory to be of the form 

^eff = J dt d^p tpP [{idt + /x)7o - p- 7 ) 'ipt (6) 

In the Wilsonian RG scaling, we eliminate all modes with energy and momenta |/co|, |^| > A, 
where A is our cutoff. As discussed above, on the remaining degrees of freedom the operator 
{idt + /i)7o - P ■ 7 0{A) and therefore scales like s. We deduce that for to remain 

invariant, ip must scale as xhe term representing the Fermi surface scales as 

indicating that it is a relevant operator and that hence the existence of the surface is natural. 
Now consider the four fermion operator 

iG J dt dVi d^p2 d^p3 d^p4 ipkX'kiipi ipiTf,ijipj (7) 

where F contains any gamma matrix and flavor or color structure. Naively, (0) is irrelevant 
since it scales as s if we assume we have I close to zero and hence 

S^{pl +P 2 -P 3 - pi) - d^{ki + k 2 - ks- k^) , (8) 

which does not scale. Higher dimension terms with extra powers of the helds or terms with 
additional derivatives are clearly irrelevant operators as well. (Again, because of the spurious 
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symmetry higher derivative terms must enter in combination with /x, and hence scale like 
s.) The only operators that survive are those satisfying the kinematic constraint pi ~ —P 2 ; 
so the delta function becomes 

<^(^1 + h — h ~ h) (9) 

and scales as (/ dl6{l) = 1). The resulting four fermion operator is marginal. 


The only interactions in the effective theory are therefore four fermion operators at a 
special kinematic point, and the model is exactly soluble by resummation. The only diagrams 
allowed by the momentum structure are the bubble diagrams found at large N in the familiar 
0(N) model. The interaction generates cooper pair formation through the gap equation (see 
Bailin and Love in |^) 

A = iG — 

Jo (27r) 

where G is the coupling, T^ any associated Dirac structure and C{k) a 4 x 4 off diagonal 
component of the 8 x 8 propagator associated with the fermion vector 


rW(p)r, 


( 10 ) 


C{p) = 


A 


( 11 ) 


i^-P7o) [A(;i-/i7o)-iA-(;i + /i7o)] 

where A = 70 AI 70 . The gap integral would diverge logarithmically as po 0 and |p|^ — pJ, 
were it not for the factors of A in the denominator of C{p) which provide an IR cutoff. 
Ignoring the Lorentz structure of the integral for the moment, we have 

d*p 


A 


-iG 


r^uv 

I A 


(27r)4 (;i - P7o) ^ (;i +/i7o) 


or, 


A 


-iGA f^uv d^k dl dko 
4 

GNA 


In 


A (27r)^ 

Ac/y' 


{ko + l - ie){ko - I + ie) 


A 


( 12 ) 

( 13 ) 

( 14 ) 


Here N = J dJk/{2TTY = /i^/27r^ is the density of states at the Fermi surface in the lowest 
order approximation. The Cooper pair condensate is of the form 


A = Af/y e G 


( 15 ) 


where c is a constant. 

Thus we obtain the result (well known in the context of non-relativistic superconductiv¬ 
ity) that any weakly attractive operator of this form generates cooper pair formation. A 
repulsive interaction [G < 0) only has the consistent solution A = 0. This conclusion in the 
relativistic case depends on the precise Lorentz matrix structure of the interaction and the 
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condensate and is the subject of study below. Although the condensate appears to be formed 
due to weak interactions this is a misleading consequence of the exact solubility of the gap 
equation. In fact as we will show in section 4 the couplings become strong as the effective 
theory approaches the Fermi surface, and from this point of view it is the strong coupling 
physics that generates the condensate. By analogy, if one could solve the Swinger-Dyson 
equations of zero density QCD at short distance, where the coupling is small, one might in¬ 
correctly conclude that the usual chiral condensate results from weak coupling effects, while 
it is in fact due to strong IR effects. 

To determine the behavior of QCD at high density we must first perform the matching 
between QCD and the four Fermi theory (section 3). We will see that several types of four 
fermion interaction are generated. In section 4 we will perform a RG analysis of how those 
couplings run to the Fermi surface to identify the dominant attractive interaction. 


3 Matching to QCD 


At very high density, /i 3> Aqcd, we expect the four quark coupling to be dominated by 
one gluon exchange due to asymptotic freedom. By taking fi sufficiently large g{fi) can be 
made arbitrarily small and one loop effects are then negligible. At tree level QCD generates 
several four fermion vertices. A pair of quarks in the 6 representation of SU{3) can scatter 
to a 6, and quarks in the 3 can scatter to the 3. However, the 6 and 3 can not mix. In 
addition because the introduction of /i breaks the 0(3,1) invariance of the theory we must 
treat the temporal and spatial components of the couplings independently. We thus have 
four interactions in the four fermi theory: 

iGl {ScaSdb + ScbSda)'ipc7o'lpa'lpd7o'lpb 


^Og i^^ca^db Scb^da}'^c'yi'^a'4^d'yi'^b 

iGl {ScaSdb - ScbSda) ^clo'tpa'tpdlo'tpb 

iG\ {ScaSdb - ScbSda) 
where we have used the decomposition 

(T^)ca{T‘^)db = -{Sca^db + ^cb^da) — - {^ca^db “ ^cb^da) 
0 3 


(16) 


(17) 


6 



Henceforth we will include signs from the contraction of spacelike 7 ^ matrices in the coupling 
constants. The one-gluon approximation yields the following matching conditions: 

/^O _ 1 Anas _ _ 1 4TTaa 

^6 ~ 6 '-^3 ~ 3 Qg'j 

_ 1 47ras _ 1 Attcis ' ' 

^6 — 6 k'2 '-^3 “ 3 fc2 ’ 

where k is the four-momentum in the gluon propagator. 

However, this is too naive and leads to non-local interactions dominated by arbitrarily soft 
gluon exchange (A: —0), with all the consequences of asymptotic freedom. It is reasonable 
to assume that the non-zero quark density screens long-range gluon exhange, leading to an 
effective gluon “mass”. At one loop order, a chemical potential leads to Debye screening of 
the temporal component of the gluon that can be calculated perturbatively |l[]. The details 
are not important here, but the result is of the form 

Ml ~ . (19) 

It is generally believed that the spatial or magnetic components of the gluons acquire masses 
only non-perturbatively. If the hnite density behavior is similar to the hnite temperature 
case, one would expect 

. ( 20 ) 

In this paper we merely assume the existence of both a magnetic and electric screening mass, 
with the electric mass somewhat larger. Our results will be qualitatively independent of the 
precise values of these masses, which provide a cutoff on infrared effects. In the four fermion 
theory we therefore introduce a form factor 

(F - Ml) Ml 1,1 + klMl) Ml 

Since /cq —>■ 0 as we approach the Fermi surface, < 0, and the largest interactions come 
from scatterings with small momentum transfer. In fact, our previous arguments concerning 
the spurious U{1)b symmetry imply that the higher order terms in this expansion are sup¬ 
pressed (correspond to irrelevant operators) near the Fermi surface. Hence in the s —0 limit 
we need only retain interactions with constant form factor (no derivatives) in our RG anal¬ 
ysis. In the BCS analysis 0 a constant potential is also used, and this is usually regarded 
as a crude approximation. However, the symmetry argument implies that it is accurate 
near the Fermi surface. Note that the symmetry argument does not imply that the 
terms can be neglected at the matching scale. Indeed, one can easily see that neglecting the 
higher order terms on the r.h.s. of (|^) leads to a large overestimate of the strength of the 
interaction. Note that we have performed the matching in Feynman gauge. In the one loop 
approximation results such as the location of the Landau pole are gauge invariant, so we 


Ml 


+ 


( 21 ) 
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can make any choice of gange. However, some gange dependence may arise from the glnon 
masses, which we have inserted by hand. 


In order to obtain an accnrate matching of a fonr fermion operator with constant form 
factor to the one-glnon interaction, we average the momentnm dependent fnnction in (^) 
over all incoming and ontgoing qnark momenta, with both qnarks near the Fermi snrface. 
This is eqnivalent to projecting ont the zero angnlar momentnm component of the operator. 
It is possible to separate any operator into eigenstates of angnlar momentnm, which do not 
mix nnder the RG flow. Since we expect that the eventnal condensate will be in the zero 
angnlar momentnm channel (see 0] for a carefnl justihcation of this in the case of the one- 
glnon exchange kernel), we do not lose any information by only retaining this component. 
Let p and q denote the incoming and ontgoing momenta of a qnark line, and 


= {p — qY = p^ + — 2p ■ q 


( 22 ) 


For p and q near the Fermi snrface, and 6 the angle between p and (f, this becomes 


= —2p^(l — cos 6*) 


( 23 ) 


Averaging the form factor (^T]) over all relative orientations of p and q, we obtain 
1 /■ o 1 1 r+i 1 




d^k 


k^-M^ 


4p2 

1 

V 


dz - 
-1 1 


z -f MV2p2 


In 


V 


( 24 ) 


We see that in this approximation the matching conditions depend only logarithmically on 
the glnon masses A more sophisticated scheme for estimating the strength of the 

constant form factor of the fonr fermion interaction might yield a slightly different relation 
than (1^) , bnt we expect it to still depend only weakly on the screening masses. The 
logarithm in is of order ln[ 5 f“^] or ln[ 5 f“^], and is somewhat larger than one. 

Onr matching conditions are then: 



The scale of this matching is of order p, and corresponds to the typical momentnm transfer 
in the scattering of two qnarks on the Fermi snrface. 























(0,p) (l%,k) (0,g) 



(0,-p) (-ko,-k) (0,-q) 


Figure 1: One loop diagram from gluon exchange which drives the RG flow of effective four 

fermion interactions. 

4 One Gluon RG Flow 

The Fermi liquid displays interesting behavior such as superconductivity because in the effec¬ 
tive theory the marginal four fermion interactions flow to strong coupling as they approach 
the Fermi surface. To display the basic behavior, consider a theory with just the simple 
interaction 

iGxjjtjj'ilj'ilj . (26) 

This interaction leads to a one loop graph which is the analog of hgure 1 with the gluon lines 
contracted. (Note that because of the restricted momentum structure of the vertices in the 
effective theory the diagram with crossed gluon lines has no counterpart and does not play a 
role in the RG flow.) It contributes a logarithmic divergence to the running in the presence 
of a Fermi surface: 

^2 f i 1 [_ i _ 

J (27r)4 -h fi'jo - + Plo - ie 

Performing the gamma matrix algebra and taking the limit ko, |/| —0, |/cp —this 
becomes 

—(7o)p'(7o)fc« + 2 (7a)p(7a)fc« I (28) 

where the log divergent part of the integral I is given by 

1 r dko dPk dl 1 

aJ (27r)"^ {ko + I — ie){ko — I + ie) 



I 
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~ -N In 
4 


A 


IR 


A 


uv 


(29) 


N = J d?k/{2'KY = /i^/27r^ is the density of states at the Fermi surface in the lowest order 
approximation. Here we have moved from an effective theory with cutoff Auv (^o, |^| < ^uv), 
to a new effective theory with cutoff A/r. As we approach the Fermi surface {Ajji —> 0), the 
coupling G runs logarithmically. 


In this fashion we can calculate the one loop beta functions for the vertices in (p! 6 D as 
a function oi t = ln(A//{/A[/v/). It is most convenient to perform the calculation (in the 
massless limit) using two component Weyl spinors. The gluon exchange conserves helicity 


and flavor quantum numbers. From hgure 1 and (^ 8 [) it can be seen that the graphs, in two 
component notation, involve the contractions of 6 sigma matrices (associated with the gauge 
vertices and resulting from the loop integral). The following a matrix identities are useful 


a' 


= 7((T*)ad(cT*),i,^ - 6(cr°)arf(cr°)a,^ 


(cr*)afe(o-^)bc(o-°)cd(cr*)x3/(cr^)2yz(cr°)^«; = -2(a*)ad(cr*)3,^„ + 3(cr‘^)ad(cr‘^) 




(31) 


(^')ab{(^)bc{(^)cd{(^')xy{(^)yz{(^)zw = {(^')ad{(^')a 


(32) 


(0-°)a6(0-°)bc(cr°)cd(0-°)xy(0-°)j/^(cr°)z«; = ad{(^^) xw (33) 

The gamma matrix algebra at one loop distinguishes the couplings between two left (right) 
handed quarks and those between left and right handed quarks. This is a result of the diver¬ 
gence in the diagram of hgure 1. Without a chemical potential this diagram is not divergent 
for the restricted shell of momenta under consideration, and this discrepancy between the 
two interactions is not seen. The RG eqns are 

y GaLL ~ ^iGaLL^) (34) 

f ( 5 ( 0211 )" - falLL GIll + f (G2ii)") (35) 

(02ir) + 2(j2ir 0.,ir — ((j.,ir) ^ (36) 

f (5(0211)"-5G2iRG;iR + i(G;iR)") . (37) 

We suppress the RR and RL equations since they are respectively identical to the LL and 
LR equations already displayed. The a subscript corresponds to either 6 or 3; note that 
again loop effects do not mix the two channels. 




aLL 


dt 




dt 


<iG 2 i„ 


dt 


dG 


aLR 


dt 
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These equations diagonalize to 


Taking 

solutions 


d{GaLL + 


aLL) 



dt 





dt 


d(ai^R 



dt 



G\lr) 


dt 


+ G'i 




^aLLJ 


= 0 


>(g: 


0 _ \2 

aLR '^aLR) 


G(0) and G* = —a;G(0) as the boundary condition at t 


(38) 

(39) 

( 40 ) 

( 41 ) 

0, we obtain the 


GlLLit) 

G^LLit) 

GaLR^'t) 

G\lr{G) 


1/ 3G(0)(l-a;) G(0)(l + 3a;) \ 

4 lvl + fG(0)(l-a;)t ^ 1 + iVG(0)(l + 3a;)t J 

1/ G(0)(l-a;) G(0)(l + 3a;) \ 

4 lvl + fG(0)(l-a;)t ~ 1 + iVG(0)(l + 3a;)t J 


1 

4 


|^G(0)(l-3a;) + 


3G(0)(l + a;) \ 

l + f^G(0)(l + a;)ty 


1 

4 


|^G(0)(l-3a;) 


G(0)(l + a;) \ 

l + f^G(0)(l + a;)t; 


(42) 

(43) 

(44) 

(45) 


We may now consider the boundary conditions from section 3. We expect a; > 1, since 
the Debye mass appears at lowest order in perturbation theory and the magnetic mass is 
a higher-order (non-perturbative) effect. The Fermi surface lies at t ——cx). In the the 3 
channel ^§(0) is positive and the LL (and RR) couplings have a Landau pole at 

t =-l/iVG3(0)(l + 3a;) , (46) 

satisfying the asymptotic relation The LR (and RL) couplings also have a 

Landau pole at t = —3/2iVG3(0)(l -f-x) but this lies closer to the Fermi surface than the LL 
coupling pole for all a; > 1. As we approach the pole in the LL (and RR) channel the LR 
(and RL) coupling becomes negligible in comparison. 

In the 6 channel G6(0) is negative and smaller than ^§(0). The interaction remains 
repulsive and the LL (and RR) couplings have a Landau pole at t = —1/3 NGq{0){x — 1). 
The LR couplings flow to constant values as the Fermi surface is approached. 

We conclude that the Cooper pairing dynamics is dominated by the attractive 3 in¬ 
teraction between two left (right) handed quarks. The condensate that forms is therefore 
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antisymmetric in color. The LL and RR channels correspond to the antisymmetric spin 0 
state and the condensate is either or . The one glnon exchange diagrams 

do not distinguish between the scalar and pseudo-scalar condensates, because the two are 
related by the U{1)a symmetry which at the level of one gluon exchange is a good symmetry. 
The dynamics is flavor blind but the flavor of the condensate is determined by Fermi statis¬ 
tics. It must be in the anti-symmetric singlet of flavor in the case of two light quarks. In 
the case of three light quarks there is the possibility of color-flavor locking . In reality we 
expect the large strange quark mass to displace its Fermi surface relative to that of the up 
and down quarks. This mismatch of Fermi momenta disfavors mixed strange - nonstrange 
condensates, except possibly at lower densities where the condensate might be quite large. 
The mismatch in Fermi momenta is roughly 

= Pf-Pfs = P- \/p^ - ml , (47) 

where p is the common chemical potential for the u,d and s quarks. This will be relevant 
when the condensates are of order or less than Api?. For p = 300 MeV and = 150 MeV, 
Ap^ ~ 40 MeV. In the absence of mixed condensates the strange quarks will condense on 
their own in the 3 channel which is antisymmetric in color. There is only a flavor symmetric 
state available in this case. Thus the LL (RR) interaction can only cause the condensation of 
the antisymmetric spin 0 and angular momentum L = 1 state which would be smaller than 
the u and d condensate. There is also potentially competition from the LR (RL) couplings 
which have a Landau pole closer to the Fermi surface and could give rise to condensation in 
the symmetric spin 1 and L = 0 state. Again this condensate, if it formed, would be smaller 
than that in the u and d sector. 


We are now in a position to estimate the size of the actual condensate, by computing the 
RG invariant scale (location of the Landau pole (^^ ) at which the 3, LL coupling diverges. 
This yields 


A = Af/y exp 


(48) 


]VG3(0) (1 + 3x)) ' 

Using the matching conditions from section 3, and the free-quark approximation to the 
density of states A, we have 


A = A[/y exp 


bvr 


1 -|- 3a; as 


— In 


4p2 


-IN 


(49) 


Note that for densities just above the expected phase boundary for chiral symmetry restora¬ 
tion, as ~ 0{1), and (|4^) predicts a relatively large condensate. As we increase the density, 
the condensate decreases. In order to compute the rate of decrease, we use the QCD beta 
function relation 


A 


QCD 


= Auv exp 


271 


bo O'siAuv) , 


(50) 
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where 6o = H ~ 2A^/3 for N flavors. Rewriting in terms of Aqcd, 


we find 


Aqcd V Ajjv / 


( 51 ) 


where 


3 _ hp 

1 + 3a; ln[|^] 

i 


Recalling that Acv ~ h; "we see that the condensate scale A decreases rapidly with increasing 
/i. Incidentally, equation (^TD again shows that the system can naturally support A ~ Aqcd, 
since at low density we expect Auv ~ Aqcd- We have checked that our results in (^91) and 
(^) are qualitatively consistent with the earlier gap equation calculations of [Q and [^. 


Finally, we note that ( ^9[) actually predicts that the condensate will eventually increase 
at asymptotically high densities, due to the logarithm in the exponent. It is not clear how 
reliable this prediction is, as it depends sensitively on the matching (averaging) procedure 
and on the precise form of the magnetic screening mass Mj. 



Figure 2: The gluonic correction to the instanton induced vertex. 

5 Instanton RG Flow 

Recently a number of groups have considered the possible role of instanton induced 
interactions in color superconductivity, particularly at low densities. In this section we 
study the combined RG flow of instanton and gluon-induced interactions. At sufficiently 
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high density, instanton effects are almost certainly negligible. As we lower the density, the 
relative size of the instanton-indnced operators at the matching scale shonld grow relative 
to the glnon operators. At snfficiently low densities we can imagine that the two types of 
interactions become ronghly comparable in strength. The RG eqnations will tell ns which 
interactions dominate the dynamics at the Fermi snrface. 

In the two-flavor case, which we analyse here, the instanton-indnced operators are of the 
form 

(ujiULdjidL - URdidRUi) ■ (52) 

We assnme that at the matching scale the instanton-indnced conpling fi:(0) is smaller than the 
glnon-indnced conpling G(0). If this hierarchy persists as t > —oo, then it is self-consistent 
to neglect loop diagrams with mnltiple instanton vertices. In this case the RG flow of k is 
determined by the effective held theory connterpart of the diagram in £g 2, which depicts 
the instanton operator renormalized by glnon exchange. In the effective theory we compnte 
the loop diagrams with one glnon and one instanton type vertex. We obtain the following 
one loop beta fnnctions, where again the snbscript a denotes either the 6 or 3 channels. 

dK 

= + . (53) 


Solving 
We hnd 


, we can directly compare the size of Ra to (the largest glnonic conpling) 
Ka(0) 4 


{1 + NG{0){1 + 3x) t)^/^ 


(54) 


^3Ll(^) ^3I,l( 0) (1 + 3a;) 

which shows that as we approach the Fermi snrface the instanton interaction becomes more 
and more snppressed relative to the glnonic operator. When the glnonic 3 interaction reaches 
its Landan pole at t = — 1/A^G(0)(1-|-3a;), the instanton interaction is completely negligible. 
This resnlt implies that as long as the instanton interaction is somewhat smaller than the 
glnon interaction at matching, the pattern of condensation will be determined solely by 
glnonic effects. 


The diagram in hgnre 2 also indnces operators with different Lorentz structnre from the 
bare vertices, as do graphs with mnltiple instanton vertices. For example, an interaction of 
the form (^ij7o7iq'L)(?K7o7i?L) is indnced. The ratio of this new conpling to behaves 
similarly to (^^ and thns vanishes as one approaches the Fermi snrface. In general, new 
vertices are always smaller than the bare conplings if we take fi;(0) somewhat less than G(0). 
However, if fi:(0) and G(0) are of the same size, as is possibly the case at lower densities, 
then the RG eqnations become considerably more complicated and the instanton vertex may 
play a role in Gooper pair formation. 
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6 Discussion 


The renormalization group provides a powerful tool for the analysis of QCD at large quark 
density. The highly constrained nature of the effective field theory describing quarks near 
the Fermi surface allows a detailed analysis of the relevant interactions. We hnd that the 
operators induced by gluon exchange lead to Cooper pair (diquark) condensation and color 
superconductivity in the 3, 0^ channels. The relevant order parameter for the two-flavor 
condensate is 

i’lai-p) C (cos(0) +fsin(0)75) ijjb{p) , (55) 

where C is the charge-conjugation matrix, (i,j) are flavor indices running from 1 to 2, and 
(a,b) are color indices. 6 = 7r/2 preserves parity, while other values break it spontaneously. 
Since the different 6 condensates are related by the U{1)a symmetry, which is good to 
all orders in perturbation theory, the gluon interactions cannot differentiate between the 
two possibilities, and it is up to interactions like instantons to determine whether parity is 
violated [^]. At high densities, where instanton effects are negligible, explicit quark masses 
also violate the U{1)a symmetry, and we have checked that the parity conserving vacuum 
is the true minimum. Nevertheless, at high density the entire family of parity violating 
and parity conserving vacua are nearly degenerate. This implies the existence of a light 
pseudo-Goldstone boson analogous to the rj', except that it corresponds to rotations in the 
diquark condensate rather than the usual quark condensate. This particle is a color singlet, 
pseudo-scalar. Due to this excitation the mass gap in the physical spectrum can be much 
smaller than the condensate scale A. 

Our calculations were performed assuming a sufficiently high density to justify perturba¬ 
tive matching to QCD. Strictly speaking, they should only be applied at densities somewhat 
larger than the critical density where the chiral phase transition is expected to occur. Nev¬ 
ertheless, they yield results which can be naively extrapolated to lower densities. Equation 
(1^ suggests that a condensate of order Aqcd is possible. 

We also studied instanton induced operators in the effective theory. We found that the 
RG ffow of such operators is slower than that of the one gluon operators, and hence that the 
latter reach their Landau pole first. This lends further confidence that the one gluon results 
can be extrapolated to lower densities. However, it is still possible that non-perturbative 
effects affect diquark condensation, depending on the relative sizes of instanton and gluon 
operators at the matching scale. 

Finally, we note that all of our calculations have been at zero temperature. A nonzero 
temperature would smear out the Fermi surface, eventually eliminating the condensates. 
Bailin and Love [@ studied the phase transition between the superconducting and normal 
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states in the one gluon exchange approximation, and found a fluctuation-induced hrst order 
phase transition. Since our results support the dominance of gluonic effects, they also suggest 
that the phase boundary is likely to be hrst order. 
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